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ABSTRACT 

A recent observation of the nearby (z=0.042) narrow-line Seyfert 1 galaxy 
RE J1034+396 on 2007 May 31 showed strong quasi-periodic oscillations (QPOs) 
in the 0.3-10 keV X-ray flux. We present phase-resolved spectroscopy of this ob- 
servation, using data obtained by the EPIC PN detector onboard XMM-Newton. 
The "low" phase spectrum, associated with the troughs in the light curve, shows 
(at > 4cr confidence level) an absorption edge at 0.86 ± 0.05 keV with an ab- 
sorption depth of 0.3 ± 0.1. Ionized oxygen edges are hallmarks of X-ray warm 



(N 



absorbers in Seyfert active galactic nuclei (AGN); the observed edge is consistent 
with H-like O VIII and implies a column density of Aoviii ~ 3x 10 18 cm -2 . 
The edge is not seen in the "high" phase spectrum associated with the crests 
in the light curve, suggesting the presence of a warm absorber in the immediate 
vicinity of the supermassive black hole which periodically obscures the continuum 
emission. If the QPO arises due to Keplerian orbital motion around the central 
black hole, the periodic appearance of the O VIII edge would imply a radius of 
~ 9.4(M/[4 x lO 6 M ])- 2 / 3 (P/[l/ir]) 2 / 3 r g for the size of the warm absorber. 

X 

Subject headings: galaxies: active - - galaxies: nuclei - - galaxies: Seyfert 
galaxies: individual (RE J1034+396) 



Introduction 



RE J1034+396 is a narrow line Seyfert 1 galaxy (NLS1; see lOsterbrock &: Poggelll985 



for definition of NLSls). It is well known for its high soft X-ray excess compared not only 



to other AGNs but a lso among the NLSls (IPuchnarewicz et al.l Il995l ; iBoller et al.l 11996 



Middleton et al.l 120071 ). While the origin of this soft excess is still not clear, the similarity of 



its soft X-ray spectrum with that of stellar black hole binaries in their 'high state' has been 
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used to postulate that RE J1034+396 harbors a comparatively low m ass black hole which is 
currently in a state of high mass accretion rate (IPounds et al.lll995l ). 



A recent X-ray observation of RE J1034+396 using the XMM-Newton satellite made 
between 2007 May 31 and 2007 June 1 showed a strong signature of quasi-periodic variability 
in the 0.3-10 keV light curve ( JGierliriski et al.l 120081 ) . The oscillations were transient in 
nature since they have never been seen in any prior observations of this source, nor did we 
see strong QPO-like variability in the 0.3-10 keV light curve obtain ed during a subsequent 
XMM-Newton observation made on 2009 May 31. It was noted by iGierliriski et al.l (120081 ) 
that the frequency of the observed oscillations was similar to what would be expected if 
the frequ encies of the quas i-periodic oscillations often seen in Galactic black hole binaries 



(see e.g. Ivan der Klid 120061 for a review on observations of QPOs in compact stellar X- 
ray binary systems) were to be scaled by the same factor as the mass ratio between the 
mass of the sup ermassive black hole in R E J1034+396 and the mass of typical stellar black 
holes. Recently iMiddleton fc Pond (120091 ) have performed detailed comparison of the timing 
properties of this XMM-Newton observation with that of stellar mass black hole binaries and 
concluded that the observed QPO in RE J1034+396 is similar to the 67 Hz QPO seen in 
GRS 1915+105 (a black hole binary system which, like RE J1034+396, also boasts a high 
L/L Edd ). This scaling of oscillation frequencies with mass would support the hypothesis that 
accretion physics scales with mass and that the supermassive black holes at the center of 
galaxies are scaled-up versions of stellar-mass black holes. 

The origin of QPOs in stellar black hole (and neutron star) binaries is not yet well 
understood. Various models have been proposed, but none can fully explain the wealth of 
observed phenomenology. QPO frequencies in stellar black hole binaries range from Hertz to 
kHz, and it is extremely difficult to carry out phase-resolved spe ctroscopy in these sy stems. 
In the sparse cases (e.g. in the source GRS 1915+105, reported by Morgan et al.ll 19971 ) where 
it is possible to "see" th e QPO in the light cur ves, the phase of the oscillations appear to 
perform a random walk. iMiller &: Homanl (120051 ) showed that for GRS 1915+105 the flux in 
the Fe Ka line (created by hard X-ray photons irradiating a relatively cold accretion disk) 
varies with the QPO oscillation phase of the 1 and 2 Hz QPOs, thus suggesting that the 1 Hz 
and 2 Hz QPOs may be linked to variable reflection. The strong variability amplitude and 
comparatively large period of the oscillations in RE J1034+396 make it an ideal source to 
carry out phase-resolved spectroscopy and probe changes in the spectral energy distribution 
between different phases. 

The origin of the soft excess in RE J1034+396 is not clear, and its time-averaged spectral 
energy distribution can be well mo deled by a wide varie ty of models including reflection from 
a partially ionized accretion disc (jCrummy et al.ll2006l ). Comptonized disc emission from a 
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low temperature disc, ioni zed partial covering, or a smeared disc wind seen in absorption (see 
e.g. lMiddleton et al.ll2009l ). Phase- resolved spectra could give important insights for breaking 
this degeneracy. Therefore one of our major goals in this work is to search for signatures of 
variable reflection (in the continuum and reflection lines) and/or variable absorption features 
(e.g. changes in the properties of O VII, O VIII edges which are the hallmark of the presence 
of warm absorber in active galactic nuclei). 



Middleton et al.l (120091 ) analyzed the energy dependence of the observed variability in the 
EPIC data and showed that the variability primarily originated in the high-energy photons. 
They also presented a spectral decomposition of the EPIC MOS spectrum, averaged over 
the entire observation. Here we analyze the EPIC PN data taken during this observation 
to test if there is any signature of spectral variations between the different phases. Due to 
limitations in photon statistics (larg ely caused by the high pileup on the CCD as described 



in §2] also see iMiddleton et al.ll2009l ) we could extract spectra with good signal-to-noise ratio 
from two phases only: the "high" phase spectrum by accumulating spectra near the crests, 
and the "low" phase spectrum by accumulating spectra near the troughs. In $2] we describe 
in detail the algorithm we used to create the phase-resolved "high" and "low" spectra. The 
spectral analysis is described in §3J and conclusions are summarized in §4J 



2. Data preparation 

The QPO in RE J1034+396 was noticed in an XMM-Newton observation starting on 
2007 May 31 (observation ID. 0506440101). The EPIC-PN detector was operated in conjunc- 
tion with the thin filter in full-frame mode during this observation. We extracted the light 
curves and spectral information of this data set from the EPIC PN detector using XMM- 
Newton Science Analysis Software (SAS; v. 9. 0.0) and following the XMM-Newton User Guide 
q Briefly, we selected data from the PN detector which satisfies the following criteria: pat- 
tern <4 (i.e. single and double pixel events only), events with pulse height between 0.2 and 
15 keV, FLAG=0 (the most stringent filter recommended to get a high-quality spectrum). 
The red line in Fig. [1] shows the 9 point moving averaged light curve extracted from EPIC 
PN data, from a circular region of 45 arcsecond radius around the source. A circular, source- 
free region, of radius 43.5 arcseconds on the same chip where the source was located, was 
chosen to extract background spectra. We used the SAS task epatplot to anlyze the event 
pattern information near the source region from the PN event file. In the absence of any 
pile up, the observed-to-model singles and doubles pattern fractions ratios (obtained from 



1 http://xmm.esac.esa.int/external/xmni_user_support/docunientation/sas_usg/USG 
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epatplot output) should both be consistent with 1.0 within statistical errors. Presence of pile 
up causes the singles ratio to decrease from unity and the doubles ratio to increase from 
unity. When a circular source region is chosen around the flux centroid of RE J1034+396, 
the observed distribution of single and double pattern fractions (as a function of energy) is 
highly discrepant from the expected model curves. The 0.5-2 keV observed-to-model frac- 
tion of singles in this case is 0.905 ± 0.003 and that of the doubles is 1.298 ± 0.007, showing 
a high pile up. Therefore we excised the PSF core where pile up is strongest. The optimal 
radius-of-exclusion was determined by slowly increasing the exclusion radius till the observed 
distribution of single and double events ma tched fairly well with th e expected curvesj. Thus 



we verified that the source region used by iMiddleton et al.l (120091 ). i.e. excluding the inner 



32 arcsec and using an annular region with outer radius of 45 arcsec, is an optimal choice. 
The 0.5-2 keV observed-to-model fraction of singles after the exclusion of inner 32 arcsec is 
0.983 ± 0.017 and the corresponding fraction of doubles is 1.096 ± 0.028. The average PN 
source count-rates in 0.3-10 keV energy range before and after excluding the PSF core were 
5.8 counts/s and 0.31 counts/s respectively. Appropriate photon redistribution matrices and 
ancillary region files were created using rmfgen and arfgen tasks. 



As noted by iGierlihski et al.l (120081 ) (also see Fig. [TJ where the 9-point moving average 
light curve computed from the 100-s binned PN data is shown by the thick red line), a visual 
inspection of the light curve clearly shows the periodicity of ~1 hour. These authors also 
noted that the variability is not strictly periodic over the entire observation, in particular 
during the initial ~25 ks of the observation. Therefore we adopted the following algorithm 
(also illustrated in Fig. [2]) to extract spectral data from the "high" (crest) and "low" (trough) 
phases. Let the 9-point moving averaged flux at time t be denoted by f(t). 

(1) In order to select a trough in the flux, we start with an initial approximation of the time 
of the trough, say t x . 

(2) We then search the interval (t{ — -P/4, t\ + -P/4), P being the QPO period, for a local 
minimum / m i n , occurring at time t min . 

(3) Next we search for the crest immediately preceding t min , occurring at ti with flux f\ Q . 

(4) Now starting from t m ; n , we search backwards in time within the interval (ti , t m j n ) for the 
time £_ when the flux reaches /_ = / min + r\ ■ (/ lo — / m i n ), with r\ = 0.4. This gives the left 
endpoint of the "good time interval" (GTI) for the trough at t min . 

(5) To find the right endpoint of the good time interval for the trough, we similarly search 
for the crest immediately succeeding t min , say tu- 



2c 



See e.g., http://www.astro.lsa.umich.edu/~dmaitra/epat.gifwhich shows the epatplot outputs for various 



exclusion radii. The numerical values of the radii are included in the file names shown on top right. Please 
note that the radii are in units of 0.05 arcsec, i.e. XMM's sky coordinates. 
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(6) As in the construction of £_, we linearly interpolate starting from £ m j n and search forward 
to find the time £ + when the flux reaches / + = f min + r\ ■ (/hi — /mm)- 

(7) We now have a good time imterval (£_,£+) for the trough of t m[n . A similar method is 
employed to estimate the GTI for the crests. 

At the boundaries of the light curve, we assumed that both the endpoints were local minima 
and computed the GTI for their respective neighboring local maxima. Looking at the light 
curve in Fig. [TJ it is clear that this allows a conservative estimation of the GTIs near the 
boundaries. 

The advantages of using this algorithm are: (1) it allows better treatment of asymmetric 
extrema, i.e. either when the maxima/minima are broad (e.g. the crest near £=0.5) or when 
the two crests (troughs) that enclose a given trough (crest) have unequal amplitudes (e.g. 
the trough near £=1.5); (2) it allows better treatment of light curves that are not strictly 
periodic but slightly wander in phase or frequency. 

The parameter r\ essentially determines the width of each GTI around any given crest or 
trough. The choice of its value of 0.4 in this case is primarily dictated by photon statistics. A 
value smaller than 0.4 would select narrower strips near the crests/troughs, and in principle 
should give a sharper distinction between spectra from crests and troughs. In reality however 
we are limited by photon number statistics and the signal-to-noise of the spectra starts 
becoming poorer for smaller r\. Increasing rj on the other hand (e.g. >0.5) causes the strips 
to overlap which is obviously not intended. Thus 0.4 is an optimal choice for r\. 

The selected high and low phases are shown in light red and light blue respectively in 
Fig. [TJ The algorithm chooses all the extrema correctly except the ones near £=1.3, 4.1 and 
6; it fails in these cases because the successive extrema occur in time scales significantly 
shorter than expected (nominally we expect two extrema to be separated by P/2). We have 
excluded these three time intervals from our analysis. The task gtibuild was used to create 
SAS-compatible GTI files for both high and low phases, and these GTI selections were passed 
to evselect task to extract the phase-resolved spectra. 



3. Analysis 



We used the ISIS (v. 1.5.0-7; see iHouck fc Denicolal 120001 ) spectroscopy package to an 



alyze the spectra presented here. In combination with additional tasks written by Mike 
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Nowalcl, ISIS offers ffexible optiongj, which made it a good choice for the analysis presented 
here. We grouped the EPIC PN data for all fits starting at >0.25 keV, with the criterion 
that the signal-to- noise ratio in each bin had to be >5. We then only considered data in 
the interval 0.25-10 keV. The photons statistics in the phase-resolved spectra, particularly 
those of the higher energy channels, are low because of (a) the low count rates after pileup 
correction, (b) keeping single and double pixel events only and using FLAG=0 to ensure 
highest quality spectrum, (c) the stringent S/N requirement in any bin, and obviously (d) 
the phase grouping. Galactic absorption was fixed to JV# = 1.31 x 10 20 cm -2 in all fits. The 
quoted errors correspond to 90% confidence limit, unless explicitly stated otherwise. 



3.1. Fits to the high phase spectrum 

A high-resolution X-ray spectrum of the heart of an active galactic nucleus, particularly 
a NLS1, has a complex shape, and physically motivated models must include contributions 



from an accretion disk (jShakura fc Sunvaevlll973l; 



ing of soft photons from the disk (Titarchukl 19941). reflection ( Ross fc Fabianlll993l ). disk 



Watarai et al. 



20001). Compton upsc atter 



winds ( Proga fc Kallmanl 12004 ; ISchurch fc Pond 120071), and also sometim es partial covering 
of the source by clumpy clouds near the source (e.g. iReeves et al.l l2008h . Furthermore, in 
comparison with most of the other active galactic nuclei RE J1034+3 96 shows a very strong 
soft X -ray excess whose origin is yet to be fully understood (see e.g., Middleton et al.l 120071 . 
20091 . for discussions). Here however, our initial question is whether the continuum itself 
changes shape between the high and low phases. Any continuum model which fits one of 
the phases well can be used to answer this question. Hence we restrict ourselves to simple 
continuum models. 

The high phase spectrum can be well modeled with a simple blackbody of temperature 
kTbb = 0.llo!o2 keV and a broken power-law with photon index Ti = 3.6 ±0.1 below 1.9 keV 
and T 2 = 1.9lo'g above 1.9 keV (also see Tabled]). Hereafter this model will be referred to as 
MhiJbh-bkn. While formally acceptable with x 2 /dof = 125.6/141, the power-law component 
in this model dominates not only the high energies, but also the low energies as shown in 
Fig. |3j Since only the peak of the blackbody contributes to the observed spectrum, replacing 



J Available at http://space.mit.edu/home/mnowak/isis_vs_xspec/download.html 



4 Such as rebinning (by minimum counts/bin or minimum signal-to-noise ratio/bin) and setting systematic 
errors on the fly, better handling of 'unfolded spectra', its support for parallel processing to take advantage 
of multi-cpu/cluster computing environment allowing fast computation of numerically extensive tasks such 
as estimation of confidence limits for model parameters. 
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the blackbody with accretion disc models does not change the fits significantly. The fact 
that the power-law dominated both l ow and high energie s when usin g a single blackbody 



was also noted in a previous work by [Pounds et al.l (119951 ) . Following iPounds et al.l (Il995h 



we tried fitting two blackbodies with different temperature and normalizations to fit the soft 
(~0.3-0.9 keV) and very soft (<0.3 keV) energies and a power-law for the high (>0.9 keV) 
energies. The best fit parameters for this model (hereafter MhiJlbbjpow) are presented in 
Table HJ and the model components and residuals are shown in Fig. [3j Note that the power- 
law model parameters T 2 for the model Mhi_bb_bkn, and T for the model Mhi_2bb_pow are 
largely constrained only by a few of the hardest energy bins, and consequently have large 
error bars. The 0.25-5 keV flux in the high phase spectrum is (1.26±0.1) x 10~ n ergs/cm 2 /s. 



3.2. Fits to the low phase spectrum 

First we test if the shape of the low phase spectrum is the same as that of the high phase. 
For this we multiply the best-fit high phase model by an overall normalization constant and 
search for minimum \ 2 allowing only this constant to change. The best-fit normalization 
constant has a value in the range of 0.90±0.02 for both models Mhi_bbJ)kn and Mhi_2bb-pow, 
when applied to the low phase spectrum. However the residuals in this case (see Fig. H|) show 
a systematic deviation in the ~0. 85-1.1 keV range, strongly suggestive of an absorption 
edge. There may also be some systematic deviation at lower energies, especially an edge-like 
absorption feature between 0.3-0.4 keV. While not statistically significant in the current 
dataset, if this indeed is the case then the continuum level at these energies is higher than 



our estimate. Some of these spectral features were also noted by iMiddleton et al.l ( 120091 ) 
in the time-averaged spectrum, but coadding spectra from both high and low phases most 
likely lowered the observability of the features (which are prominent only during the low 
phase). 

Noting that the 0.8-0.9 keV feature could be due to an H-like o xygen VIII K e dge (at 



rest frame energy of 0.87 keV, seen in many Seyfert 1 sources (see e.g. lReynoldslll997l ) due to 
presence of ionized gas in the line of sight), we next attempt to test if this feature (present 
only during the low phase) could indeed be due to the oxygen absorption edge. For this we 
included a multiplicative edge component to our models. The best-fit values after including 
an edge in both Mhijbbjbkn and MhiJ2bbjpow are given in Table [Q 

The improvement in x 2 fit statistic for model MlojbbJ)kn_edge between r = (no edge) 
and r = 0.29 (best fit) is A% 2 = 19.4, implying a single parameter confidence limit of 4.4a. 
The improvement in \ 2 n ^ statistic for Mlo.2bbjpow-.edge between r = (no edge) and 
r = 0.23 (best fit) is A% 2 = 16.6, implying a single parameter confidence limit of 4.07cr. 



Thus both low phase models need an edge at >4<r level, strongly suggesting that the 
edge is present in the data. While fits to both models MloJbbJbkn and Mlo_2bb_pow are 
formally acceptable, there remain hints of some systematic trend in the residuals shown in 
Fig. H]at lower energies. This is discussed in greater detail in §U The 0.25-5 keV flux in the 
low phase spectrum is (1.14 ± 0.1) x 10~ n ergs/cm 2 /s. 

Next we attempted to estimate the strength of the same OVIII absorption edge in the 
high phase data, by adding an edge to both the models MhiJ>b_bkn and MhiJ2bb_pow and 
redoing the fits. We found that inclusion of the edge did not improve the statistical quality of 
the fits. The OVIII edge could not be detected (at > 3<r level) in the data using either of the 
continuum models, and this confirms that the optical depth of the edge varied significantly 
between the high and low phases. 

Next we modeled the edge in the low phase s omewhat more physically using the XSTAR 



photoionization code ( iKallman fc Bautistall200ll ). A grid of XSTAR models was created from 



the following variables: (1) the column density (Nh, allowed range 10 20 - 10 23 cm -2 ), (2) 
ionization parameter (log(t;) = log(L/(nR 2 )), allowed range 2-4), (3) oxygen abundance 
(relative to solar, allowed range 0.1-4), (4) neon abundance (relative to solar, allowed range 
0.1-1). A multiplicative table model was created from this multidimensional grid. Turbulent 
velocities of 300 km/s and above cause absorption lines to become more prominent, and lower 
turbulent velocities are required to fit the data well. The best-fit parameters to the low 
phase data using this model in conjunction with a continuum model composed of double- 
blackbody+power law suggest iV#/(10 21 cm -2 ) = 4 ± 1, log(£)=3.2+Q; 2 , and a relatively 
high oxygen abundance (2.11°;!). A sub-solar neon abundance (0.13) is preferred by the 
best-fit solution but its error-bar could not be constrained. See Table [1] for the detailed fit 
parameters and statistics. The best-fit continuum parameters in this case were very similar 
to those obtained using the single edge model (MloJZbbjpow). 

It is possible that oscillations in the flux of the source cause the ionization of intervening 
gas to vary, and that the gas is too ionized to give significant OVIII absorption during high 
phases. This possibility can be tested using our XSTAR model for the warm absorber gas. If 
this were the case, then it ought to be possible to model the high phase spectrum using the 
same warm absorber abundances as obtained for the best-fit to the low phase data, but with 
a higher ionization fraction (~10-20% higher, which is typically the flux difference between 
the crests and peaks). However the high phase spectrum cannot be modeled with 10-20% 
enhanced ionization fraction. In fact, detailed modeling shows that reasonable fits to the 
high phase spectrum can be obtained only if the ionization fraction is a factor of three higher 
than that required by the low phase spectrum. Thus it seems unlikely that the changes in 
the high and low phase spectrum are driven solely by flux viariations in the source. 
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We have also repeated the entire analyses using a different, somewhat larger exclusion 
radius of 35 arcseconds (smaller exclusion radii were not considered as they would be affected 
by pile up). In this case the errors are larger due to poorer photon statistics, but the 
conclusions are same, i.e. an edge is detected (at > 3a) level in the low phase, but no edge 
is detected in the high phase. Similarly we repeated the analyses extracting only pattern=0 
events (single pixel events), and excising the PSF core of radius 30 arcseconds. In this case 
also the OVIII edge was detected at 4.1<r confidence level in the low phase spectrum whereas 
no evidence of the OVIII edge could be found the high phase spectrum. Thus our main 
observation, viz. that the OVIII absorption edge is only seen during the troughs in the 
oscillations, and not during the crests, appears to be quite robust and is not some form of 
detector artifact. 



Discussion and Conclusions 



In this work we have presented phase-resolved spectroscopy of the narrow line Seyfert 1 
galaxy RE J1034+396. We have used the XMM-Newton observation of this source between 
2007 May 31 and 2007 June 1, when the X-ray light curve showed strong quasi-periodic oscil- 
lations. We have extracted spectra during the high (crest) and low (trough) phases. Simple 
continuum models assuming either a blackbody plus broken power-law, or two blackbodies 
plus a power-law model, fit the high phase spectrum quite well. The required low-energy 
power-law index is very steep (r ~ 3.6 ± 0.1) for a single blackbody plus broke n power-law 



model, which is typical for narrow line AGNs w ith a high soft X-ray excess (jBoller et al. 



1996; Pounds et al. 



19951 ; iMiddleton et al.l 120071 ). The spectrum can also be well modeled 



using two blackbodies and a power-law continuum. Due to the steep decline of the spec- 
trum, and high pileup on the PN detector which severely limits the photon statistics, we 
do not have good quality spectral data beyond ~2-3 keV. Therefore for these spectra we 
cannot very well constrain the relatively flatter 2- 10 keV continuum usually seen in this 
source (jCrummy et al.ll2006l ; IMiddleton et al.ll2009i ). although fits to a two-blackbody plus 
power-law model show that the > 2 keV spectrum is comparatively flatter with a photon 
index of ~ 2. 

The low phase spectrum cannot be adequately described by scaling the overall normal- 
ization of either of the best fit high phase models, suggeting that the shape of the low phase 
spectru m is different from that o f the high phase. Analyzing the r.m.s. variability at different 
energies IMiddleton et al.l (120091 ) also concluded that the variability is mainly from the hard 



photons. The right panel of Fig. H] shows that the differences between the low and high phases 
are most prominent at higher energies. In particular, the residuals in Fig. H] show a sharp 



-10- 



drop in the continuum flux between 0.8-0.9 keV, which then slowly increases toward the ex- 
pected continuum flux at higher energies, reminiscent of an absorption edge. The inclusion 
of an edge improves the fits significantly (at >4<r level) for both single-blackbody+broken 
power-law model and double-blackbody+power-law model. The best fit values for the edge 
threshold energy and absorption depth are 0.87 ± 0.04 keV and 0.29 ± 0.11 respectively for 
the single-blackbody+broken power-law model, and 0. 861^06 keV and 0.23io!ii respectively 
for the double-blackbody+power-law model. The edge is most likely associated to the 0.87 
keV (rest frame) K-edge of (H-like) oxygen VIII which is quite commonly seen in Seyferts 
(see e.g., lReynoldslll997l ). and is indicative of the presence of optically thin, photoionized 



matter along the line of sight. RE J1034+396 however is at a redshift of 0.042; therefore the 
0.87 keV edge (in rest frame) should appear at 0.84 keV (easily within the 90% confidence 
error bar for both models). The presence of a warm absorber should create other features 
in the spectrum, e.g. the 0.74 keV K edge from He-like O VII. While the residuals in Fig. 
do show a drop in the flux, it is not very significant statistically. 

To estimate the column d e nsity (iVovm) of the O VIII ions responsible for the edge 
we used Eq.(l) of IVerner et al.l (119961 ) and calculated the absorption cross section a (for O 
VIII, a[E = 0.87 keV] = 0.098 Mb). Since r = JV vm • <r, we obtain JVovm = 3 x 10 18 
cm -2 . Estimating the number density (^ovin) requires a size-scale for the warm absorber. 
Assuming that the QPO originates in the accretion disk, the warm absorber is could be 
confined within the orbit of the QPO, or in case of an outflow it originates inside the QPO 
orbit . The mass of the supermassive black h ole in RE J1034+396 is in the range of 1 — 7 x 10 6 
M Q ( Bian &: Huang||2010l ; IZhou et al.ll2010l ). A circular Keplerian orbit would have a radius 



of 9.4(M/[4 x 10 6 M o ])~ 2//3 (P/[l/ir]) 2//3 r g , which would be an upper limit on the size of the 
absorber in this case. This would give novm = 3 x 10 5 cm -3 for the number density of OVIII 
ions. Note however that this assumes a Keplerian orbit origin for the QPO, and therefore 
depends on the mass of the central supermassive black hole. Thus in RE J1034+396 we may 
be seeing an infalling "blob", the emission from which is periodically absorbed as it passes 
behind a warm absorber in the immediate vicinity of the central engine. If the "blob" is at 
~9 r g , that would imply the ionizing flux must be generated within 9 r g (or else we would 
not see absorption), and this would be easier to accomplish if the black hole is spinning since 
the disk can get down to a s maller radius (e. g . as l ow as 1.25 r g for a maximally spinning 
black hole) . Recent work by ICrummy et al.l (120061 ) has shown that NL S1 spectra c an be 
modeled using reflection from a disk around a spinning black holes (see iMillerl 120071 . for a 
review) . 

We also created a multidimensional grid of models using the XSTAR photoionization 
code, to create a more physical model for the warm absorber and explore a wide range 
of column density, ionization parameter, oxygen and neon abundance. Fits to the low- 
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phase data using this model suggests a warm absorber column density of ~ 4 x 10 21 cm -2 , 
log(£)=3.2+Q; 2 , and an oxygen abundance of 2.1+q;® relative to solar. The neon abundance 
could not be well constrained from the data, although the best-fit solution suggests a sub- 
solar value. 

Another possible physical scenario for the difference between high and low phase spec- 
trum could be that the warm absorber is along the line of sight, between us and the source 
(where the QPO originates), and the warm absorber responds to flux variations in the source 
itself, i.e. it is more ionized and transparent at high fluxes and less at lower fluxes. However, 
based on our modeling of the spectra using physically motivated models generated using 
XSTAR this scenario appears unlikely because the change in luminosity of the ionizing flux 
between low and high phases is too small to create the observed difference in the depth of 
the OVIII edge (also see §3J. 

It is interesting to note in this context that a som ewhat complementar y situation has 
been observed in the Seyfert galaxy NGC 1365, where iRisaliti et al.l (120091 ) have reported 
a possible transit of an obscuring cloud (with JV# ~ 3.5 x 10 23 cm -2 and other inferred 
properties similar to that of a broad-line region cloud) in front of the central X-ray source. 



Signatures of variable absorption on short timescales have been previously observed in 



high resolution Chandra spectra of the stellar black hole candidate H1743-22 (IMiller et al. 



20061 ). For H1743-22, the lines were slightly blueshifted, suggesting an outflow origin. For 
RE J1034+396, we do not have definitive evidence for an outflow, but the best fit edge energy 
for both models is higher (though within 90% confidence interval) than the rest frame energy 
of the oxygen edge. If the warm absorber in RE J1034+396 is indeed in an outflow, then the 
observed oscillations could originate due to an instability in the inner acc retion disk threaded 
by a poloidal magnetic field , giving rise to an outflowing wind or jet ( JBlandford fc Payne 
19821 ; iTagger fc Pellatl 119991 ). The observed oscilla tions in this case would correspond to 
scaled up versions of low- frequency QPOs (see e.g.. Ivan der Klia l2Q06h seen in stellar mass 
X-ray binaries. 

Since the presence of a warm absorber is bound to have noticeable imprints on other re- 
gions of the electromagnetic spectrum, especially in UV and soft X-rays, future deep UV and 
X-ray observations of RE J1034+396 will be useful not only in constraining the origin of the 
oscillations, but also in gaining a better understanding of the physical emission mechanism 
of sources like RE J1034+396 which show a strong soft-X-ray excess. Future X-ray missions 
with i mproved sensitivity and larger collect ing area like the International X-ray Observatory 
(IXO; IMiller et al.ll2009t IWhite et al.ll2010l ) would play key role in understanding accretion 
geometry as well as emission/absorption mechanism in these sources. 
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Table 1. Summary of fits. 



Phase 



Model 



Best-fit parameters 2 Best-fit statistics 



High 



MhiJbbJbkn 



kT _ n i +o.o3 

lil bb — u - 1 -0.02 

JV 66 /10~ 5 = 3.2^ 
Ti = 3.6 ±0.1 
Ar p /10- 3 = 1.5 ±0.1 



h2.9 
-1.7 



E b 

r 2 



1.9 

1.9 



+0.6 

-0.3 

+0.4 

-0.9 



High Mhi_2bb_pow kT bbl = 0.06 ± 0.01 

AW/10- 4 = 2.8+Q; 4 
kT bb2 = 0.14 ±0.02 
N bb2 /W- 5 = 6.9tH 
T = 2.0 ±0.4 
N p /W- 4 = 6±1 



X 2 /y = 125.6/141 



x 2 /u = 130.0/141 



Low 



MloJbb-bkri-edge kT\ 



bb 



0.09 



-0.01 
-0.01 



2.9 



+0.5 
-0.4 



N bb /10- 5 
N p /W- 3 = 1.3 ±0.2 
Ti = 3.6 ± 0.03 
E b = l-8_ ; 2 
T 2 = 2.3 ±0.5 
E e = 0.87 ±0.04 
t = 0.29 ±0.11 



x 2 /u = 136.0/125 



Low 



Mlo_2bb-pow-edge 



06 +0 - 01 
u - UD -0.02 



2 4 +cu 
z - 4 -o.o 



N bbl /10- 
kT bb2 = 0.14 ±0.04 
N bb2 /W~ 5 = 5.6 ± 0.2 
T = 2.36" 1 



N p /W 
E e = 0.86 
r = 0.23 



HX06 

-0.58 

6.7 

+0.05 
0.06 
0.12 
0.11 



+0.4 
-2.7 



X 2 /y = 133.6/12* 



Low MloJlbbjpowjxsgrid kT bb \ = 0.06 ± 0.01 
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Table 1 — Continued 



Phase Model 1 Best-fit parameters 2 Best-fit statistics 



JWio- 


-4 _ 9 4 +0.1 

— z - 4 -o.o 


kT bb2 = I 


114 ±0.01 


N bb2 /10- 


5 = 5.6+ - 2 


r = 2.42 


+0.01 
-0.51 


iv p /io- 4 


= 7.4 ±0.4 


N H /W 21 


= 4±1 


log(C) = 


q 9+0.2 
°- z -0.1 


OAbund = 


_ 9 i +0.8 
" Z - ± -0.5 



Ne A bund = 0.14 



+0.86 
-0.04 



X 2 /v = 134.3/126 



1 In ISIS/XSPEC notation the model definitions 
are: Mhi-bbJbkn = phabs (zbbody + bknpower) ; 

Mhi.2bb.pow = phabs (zbbody (1) + zbbody (2) + 

powerlaw) ; MloJbb_bkn_edge = constant [phabs 

(edge (zbbody + bknpower))]; MloJ2bb_pow.edge = 

constant [phabs (edge (zbbody(l) + zbbody(2) + 
powerlaw))]; MloJlbbjpowjxsgrid = constant [phabs 
(xsgrid (zbbody (1) + zbbody(2) + powerlaw))] 

2 The parameter kT^ is the blackbody temperature in keV. 
The blackbody normalization N b b = L3g/[.Dio(l + z)} 2 , where 
L39 is the source luminosity in units of 10 39 erg/s, z = 0.042 
is the redshift, and D\q is the distance to the source in units 
of 10 kpc. N p is the number of power-law photons/cm 2 /s/keV 
at 1 keV, and T is the photon index of the power-law. When 
using a broken power law model Ti is the photon index for 
E < E b , and T 2 is the photon index for E > E b , where E b is 
the break energy in keV. The absorption edge is parametrized 
by its threshold energy (E e , in keV) and absorption depth at 
the threshold (r). The model xsgrid is a multiplicative table 
model created from and XSTAR grid with following variables: 
column density (Njj), ionization parameter (log(£)), oxygen 
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abundance relative to solar (OAbund)> and neon abumdance 
relative to solar (NeAbund)- See $3] for details of the model. 
Galactic absorption was fixed to Njj = 1.31 x 10 20 cm -2 in all 
fits. The errors correspond to 90% confidence limit. 
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4 5 

Time (10 ks) 

Fig. 1.— The red line shows the 9 point moving averaged light curve extracted from 
EPIC PN data, from a circular region of 45 arcsecond radius around the source, i.e. before 
excising the central, piled up core of the PSF. The start time here corresponds to 2007 May 
31, 20:11:01 UTC. The high and low phases selected using the algorithm described in §2] are 
shown in light red and light blue respectively. Please see the electronic edition of the journal 
for a color version of this figure. 




[1] Initial guess for t min 

W i'miir 'm\w 

[3] (Ho, f| ) 

[4](t,fJ;f.-f m , n + 0.4(f lo .f min ) 

[5] (t hi , W 

[6](t + ,f + );f + = f min + 0.4(f hi -f min ; 

[7] (t_, t + ) : GTI for the trough 



Fig. 2. — An example showing the selection of good time interval near a trough, following 
the algorithm described in §2j 
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Fig. 3— Left panel: Best fit blackbody+broken power-law model (Mhi_bbJ)kn in Tabled]) 
and residuals to the high phase spectrum of RE J1034+396. The power-law component is 
shown by the blue histogram, the blackbody in green, and the total in red. Right panel: 
Best fit two blackbodies and a power-law model (MhiJlbbjpow in Table [I]) and residuals 
to the same dataset shown in the left panel. The power-law component is shown by the 
blue histogram, the high temperature blackbody in green, the low temperature blackbody 
in violet, and the total in red. Please see the electronic edition of the journal for a color 
version of this figure. 





^ O : 



Energy (keV) 



Energy fkeV) 



Fig. 4.— - Left panel: Best fit to the low phase data using the model MhiJobJokn multiplied 
by an overall normalization constant and allowing only the normalization constant to vary, 
thus keeping the spectral shape same as MhiJabJokn. Right panel: Same as left panel but 
for model MhiJlbbjpow. In both cases note the systematic deviations in the residuals in the 
~0. 85-1.1 keV range. 
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Fig. 5. — Best fit to the low phase data using models MloJ)bJ)kn-edge (left panel) and 
Mlo-2bbjpow-edge (right panel) which includes the (H-like) oxygen VIII K edge. The best 
fit parameters are given in Table [TJ 



